Distributions of Waiting Times of Dynamic Single-Electron Emitters 
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The distribution of waiting times between elementary tunneling events is of fundamental impor- 
tance for understanding the stochastic charge transfer processes in nanoscale conductors. Here we 
investigate the waiting time distributions (WTDs) of periodically driven single-electron emitters and 
evaluate them for the specific example of a mesoscopic capacitor. We show that the WTDs provide 
a particularly informative characterization of periodically driven devices and we demonstrate how 
the WTDs allow us to reconstruct the full counting statistics (FCS) of charges that have been trans- 
ferred after a large number of periods. We find that the WTDs are capable of describing short-time 
physics and correlations which are not accessible via the FCS alone. 
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Introduction. — Investigating the electrical fluctuations 
in a nanoscale conductor is an attractive method to probe 
and characterize the physical mechanisms and correla- 
tions that determine a given quantum transport process 
[l[ . In one approach, the stochastic number of transferred 
particles during a long time interval, the so-called full 
counting statistics (FCS), is studied [3. FCS already has 
a significant history in the theory of mesoscopic physics, 
and recent measurements of current fluctuations in sub- 
micron devices have shown that FCS is no longer just an 
interesting theoretical concept Q- It is also becoming 
an important experimental tool to examine interaction 
and coherence effects in nanoscale systems under out-of- 
equilibrium conditions. 

The FCS, however, does not provide a complete picture 
of a charge transport process on all relevant time scales. 
Conventionally, FCS is defined in the long-time limit, 
where a large number of charges have passed through 
the conductor. Only the zero-frequency components of 
the current fluctuations can then be addressed and im- 
portant short-time physics may be lost. Recently, sys- 
tematic theories of finite-frequency FCS have been de- 
veloped in order to calculate frequency-dependent noise 
and higher-order cumulants but at this point these 
methods are still restricted to systems without any ex- 
plicit time dependence. An alternative and particularly 
intuitive characterization of the charge transfer process is 
provided by the distribution of time delays between sub- 
sequent physical events, also known as the waiting time 
distribution (WTD). 

While WTDs have been studied intensively in other 
fields of science, e. g., in single molecule chemistry [B, HJ 
and in quantum optics 0, [q, they have only received 
very limited attention within the field of quantum trans- 
port. Exceptions include a few theoretical works on non- 
driven systems > but a conceptually simple example 
of a nanoscale electronic system where the usefulness of 
WTDs is clearly demonstrated has to date been missing. 
In this Letter, we show that the WTDs of periodically 
driven single-electron emitters, Fig. [IJ provide a very 
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FIG. 1: (Color online) Single-electron emitter and waiting 
times, (a) Nanoscale system coupled to source (S) and drain 
(D) electrodes. Only a single electron at a time can occupy 
the nanoscale system, e. g., because of strong Coulomb in- 
teractions. Unidirectional transport takes place from source 
to drain due to the periodically modulated rates Ts(i) and 
r_o(t). (b) Current pulses in the source (blue [dark gray]) and 
drain (red [light gray]) electrodes together with illustrations 
of the corresponding waiting times At (dashed lines) between 
elementary tunneling events, absorption (a) and emission (e). 



useful view on the charge transfer statistics and corre- 
lations in such systems. In particular, we evaluate the 



WTDs of a mesoscopic capacitor [12h14| | , which serves as 
a prime example of the usefulness of WTDs. We derive 
expressions for the WTDs which are applicable also to 
certain types of quantu m p umps IB" 12 1 and nanoelec- 
tromechanical systems [18j. We demonstrate how the 
WTDs not only allow us to reconstruct and obtain the 
FCS of emitted electrons, but additionally they contain 
information about the charge transfer process on short 
time-scales which is not available in the FCS alone. As 
we show the WTDs describe the charge transfer process 
on all important time scales. 

System. — We focus on systems consisting of a source 
and a drain electrode coupled to a nanoscale conductor, 
Fig. [TJa), biased such that single-electron transport is 
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unidirectional from the source to the drain. The tunnel- 
ing rates to and from the conductor, r,s(t) and To{t), 
respectively, are time dependent. The probability P(t) 
for the conductor to be occupied by an electron obeys 
the evolution equation 

d t P(t) = -T D (t)P(t) + T s (t)[l - P(t)] , (I) 

where 1 — P(t) is the probability for the conductor to 
be empty. This model suffices to illustrate the basic con- 
cepts of WTDs which are of interest here. We concentrate 
on systems, similar to the ones mentioned in the Intro- 
duction, where the tunneling rates are periodic in time, 
such that T a (t) = T a (t + T), a = S,D, with T being the 
period. 

Waiting time distributions. — We consider the wait- 
ing times between different tunneling events, Fig. [TJb). 
These consist of events, where the conductor either ab- 
sorbs an electron from the source or emits an electron 
via the drain. Because of the probabilistic nature of 
the charge transfer, the waiting time At between such 
events is itself a random variable. We use u; ea (At,t ) 
[w ae (At, to)} to denote the probability for the first emis- 
sion event to occur at time At + to given that absorp- 
tion occurred at the random earlier time to [and simi- 
lar for absorption following emission]. The same defi- 
nitions apply to the WTDs for tunneling events of the 
same type, w ee (At,to) and w aa (At,to). Since the oc- 
cupation of the conductor is either or I, two events 
of the same kind cannot happen simultaneously and 
w ee (0,to) = if 00 (0,to) = for all to- For nondriven 
systems, where the tunneling rates are time-independent, 
translational invariance with respect to time implies that 
the WTDs do not depend on t 0. 

To calculate the WTDs we first express the source and 



drain mean currents at time to as (Is (to)) — ^s(to)[l — 
P(t )} and (I D (t )) = T D (t )P(t ). The currents are 
proportional to the probabilities of absorbing and emit- 
ting an electron, respectively. Additionally, we need 
the conditional currents, e. g., (/£, (At, to)) = To (to + 
At)P a (At,t ). Here, P a {At,t ) is the survival probabil- 
ity of an electron at time to + At given that it was ab- 
sorbed at time to, such that P a (0,to) = 1. The survival 
probability F a (At,t ) obeys Eq. (U with T s (t) = 0, i. 
e., d At P a (At, to) = -T D (t + At)P a (At, t ). The WTD 
between absorption and emission is now w ea (At, to) — 
N ea (Is (to)) (I'd (At, to)), where Af ea is a normalization 
constant. The WTD is T-periodic in the second time 
argument, w ea (At, to + T) = w ea (At,to). This allows 
us to consider only the finite time interval to € [0, T] 
and choose the normalization constant Af ea , such that 
J T dto Jp°° d(At)w ea {At, t ) = 1. 

At this point, the WTD depends not only on the wait- 
ing time At, but also on the absolute time to at which 
absorption occurred. We are only interested in the wait- 
ing time itself and therefore integrate out to- Defining 
f = J dto f (to) for a T-periodic function /(to), we find 
for the WTD between absorption and emission 

w ea (At) =A4a(/s(to)></£(At,t )>. (2) 

Similar reasoning for the WTD between emission and 
absorption leads to the equivalent expression 

w ae (At) = Af ae (In(t )}(I e s (At,t )), (3) 

with the conditional current (T|(At,t )) = Tg(t + 
At)P e (At,t Q ). Here P e (At,t ) is the survival probabil- 
ity of the empty state. Proceeding along the same lines 
for the WTDs between events of the same kind, we find 



«w (At) =Af aa / d{At')(I s {to))(Ih(At',to))(I s (At - At', t + At')), 
Jo 

w ee (At)=M ee / d(At')(I D (t ))(I%(At',t Q ))(I D (At- At',t + At')), 
Jo 



(4) 



where the constants J\f aa and J\f ee ensure the normaliza- 
tions J °° d(At)W aa {At) = 1 and J °° d(At)uJ ee (At) = 1, 
respectively. In the expression for w aa ( At) \w ee (At)], At' 
is the time interval between the first absorption [emis- 
sion] event and the intermediate emission [absorption] 
event which finally is followed by the second absorption 
[emission] event after the total waiting time At. 

Mesoscopic capacitor. — We now illustrate our findings 
in terms of a specific example: a mesoscopic capacitor 
consisting of a nanoscale cavity coupled to external reser- 
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voirs via a quantum point contact When subject 

to fast periodic gate voltage modulations, the capacitor 
can absorb and re-emit single electrons at giga-hertz fre- 
quencies. The system can be described by Eq. (JTJ) taking 
Ts(t) = r and r^(t) = in the first half of the period, 
and T s (t) = and T D (t) =T in the second half [liGl]. 
The tunneling rate T can be controlled experimentally 
by adjusting electrostatically the opening of the quan- 
tum point contact. 

Following the steps described above we obtain a simple 
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FIG. 2: (Color online). Waiting time distribution (WTD) for the mesoscopic capacitor. We show the WTD between subsequent 
emission events W ee (At) for different values of the tunneling rate F in units of the inverse period of the driving T -1 . The mean 
charge emitted per period (the mean current) ((/)) is also indicated. The analytic result, given by Eq. ©, is compared to 
numerical simulations of the charge transport. For large values of F, (a), the charge transport is highly regular and periodic 
with the mean waiting time equal to the period, {(At)) — T. As T is reduced, (b), several peaks appear in the WTD. For even 
smaller values of V, (c) and (d), the mean waiting time is much larger than the period, ((At)) 2> T, and the overall shape of the 
WTD is determined by the envelope curve shown in blue. The mean charge emitted per period is then much smaller than 1. 



expression for the WTD between emission and absorption 
events 

W ae (At) = -^—y {e~\M - e 2 el^l} , (5) 

where e = e~ rT / 2 , f oe = r[Af - ({At/T\ + 1/2)T], and 
[x\ is the integer part of x. The result contains two 
independent structures: an internal structure (in curly 
brackets) which is periodic with T and an envelope (given 
by e L At / T J ) which is responsible for an exponential decay 
of the WTD for large waiting times At. For the WTD 
between two emission events we find 
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e 2 e^ 



'} (6) 



with £ ee = r[At - (LAt/TJ + 1)T]. Again, the WTD 
consists of an envelope function and an internal structure. 

Our analytic results are confirmed by numerical simu- 
lations of the mesoscopic capacitor, Fig. [5] For large tun- 
neling rates, Fig. E^a), the transport process is predom- 
inantly regular and periodic with one electron emitted 
almost every cycle. The WTD has a single peak centered 
around the period At ~ T. The peak, however, is not 
sharp due to the jitter in the emission process, causing 
phase noise [14j . As the tunneling rate is decreased, Fig. 
E5b), a more complicated structure appears with several 
equidistant peaks separated by the period. Two emis- 
sion events must be separated by at least half a period 
implying that w ee (At) = for At < T/2. Interestingly, 
this short-time behavior is not visible in the noise spec- 
trum of the capacitor found in Ref . [l9| . For even smaller 
tunneling rates, Figs. 0c) and[2jd), the charge trans- 
port becomes increasingly random, and subsequent elec- 
tron emissions are typically separated by several periods. 
The current fluctuations are then shot-noise like and the 
overall shape of the WTD is determined by the envelope 



function of the approximate form Ate~ < - r / 2 - |At . This cor- 
responds to the WTD for the case, where both tunneling 
rates are kept constant as r§(t) = rjj(t) = T/2 for all 
t. As we show below, much of this information is not 
available in the FCS alone. 

Full counting statistics. — To connect the WTDs to 
the FCS, the probability V(n, N) of emitting n electrons 
during a large number of periods N, we assume that max- 
imally a single electron can be emitted during a period. 
This is the case for the mesoscopic capacitor considered 
above. We can then write the probability distribution as 
V(n,N) = Y, mi ,...,rn rl ^^{rni)---w ee (rn n )5 mi+ ... +mn ^ Nl 
where w ee (m) is the coarse-grained WTD for the num- 
ber of periods m between two subsequent emission events. 
The Kronecker delta 5 m ,N expresses the constraint that 
the sum of periods between emission events m\ + . . . + m n 
must equal the total number of periods N. We have as- 
sumed that the counting of emitted electrons starts right 
after an emission event, but the specific choice of ini- 
tial condition is not important for the FCS after a large 
number of periods. 

Next, we introduce the cumulant generating function 
(CGF) S(x,N) = logX;„^(«,^)e lxn whose derivatives 
with respect to the counting field \ at x — yield the 
cumulants of n as ((n k )) = d k x S(x, N)\ x ^ . Addition- 
ally, we define the discrete Laplace transform of w ee (m), 
w ee {z) — J2 m ^ee( m )e~ mz and the corresponding CGF 
of m, Wee{z) = log w ee (—2), which similarly delivers the 
cumulants of m as ((m k )) = d k W ee (z)\ z ^ . We can then 
express the CGF of n as 
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dz 



1 _ e iX + Wee(-2) 



(7) 



In the large- N limit, the integral is determined by the 
particular pole z = Zq(x) with Zq(0) = that solves the 
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FIG. 3: (Color online). FCS of the mesoscopic capacitor. The 
normalized current cumulants Fk = {{I k )) / ({I}} , k = 2,3,4, 
were obtained from numerical simulations (FCS) as well as 
from Eq. © using both numerical (WTD) and analytical 
(Analytics) results for the coarse-grained WTD. Dashed lines 
indicate values of TT corresponding to the panels of Fig. [5] 



equation 



iX + W ee (-z) =0, 



(8) 



such that S(x,N) — > Nzq{x)- The electron current 
/ = n/N is the number of emitted electrons n over the 
number of periods N, and Zo(x) thus generates the cu- 
mulants of the current, i. e., {{I k )) = d\ ' Zq(x)\ x -+Q- For 
the mesoscopic capacitor, the large- N limit is reached for 
N » max{l,l/rr}. 

Equation ((SJ demonstrates an important and intimate 
connection between fluctuations in the current of emitted 
electrons I and the number of periods to between emis- 
sion events. In general, the equation may be difficult to 
solve for z = zq(x), but it provides us with a simple and 
systematic way of relating the current cumulants to the 
cumulants of to: taking consecutive derivatives of the left 
hand side with respect to the counting field x evaluated 
at x — 0, we find for the average current ((/)) = l/((m)) 
and the (normalized) current cumulants Fk = ((I k ))/((I)) 
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For non-driven systems, the number of periods m should 
be replaced by the continuous waiting time At. We 
then recover the known relations {{IJj = 1/ {{At)) and 
F 2 = {{At 2 )) I {{At)) 2 d, see also Ref. [6:]. Our derivation, 
however, allows us to determine current cumulants of any 
order. 

In Fig. [3] we show the FCS for the mesoscopic capaci- 
tor. We performed separate numerical simulations of the 



FCS and the WTD, and for comparison we then used 
Eq. ([9]) to obtain the normalized current cumulants Fk 
from the coarse-grained WTD. Additionally, from Eq. 
(|6|) we found analytically the CGF of to, W ee {z) = 
z + 21og[(l - e)/(l - ee z )}, and again used Eq. © to 
obtain the Fl-'s. The figure confirms the validity of Eq. 
([9]) and clearly illustrates that the FCS can be obtained 
from the coarse-grained WTD. Importantly, the proce- 
dure cannot be reversed: The WTD cannot be obtained 
from the FCS. Moreover, comparing Figs. [2] and [3l sub- 
stantial information about the charge transfer process is 
obviously lost in the FCS. In the phase noise regime, 
TT = 10 [corresponding to Fig. Ufa)], the current cu- 
mulants are close to zero due to the regular emission of 
electrons. However, contrary to the WTD, the cumulants 
are not sensitive to the jitter in the emission process. As 
the tunneling rate is lowered, FT — 2 [Fig. HJb)], several 
peaks appear in the WTD, but this is also not visible in 
the FCS, neither is the fact that two emission events must 
be separated by at least half a period. In the shot-noise 
regime, TT = 1,0.2 [Figs, [^c) and^d)], the cumulants 
approach the limiting values Fk — > (l/2) fe_1 correspond- 
ing to a Poisson process with an effective charge of 1/2. 
This is also a very different characterization compared to 
the one provided by the WTDs. 

Conclusions. — We have shown that the distribution 
of waiting times between elementary tunneling events is 
a useful tool to probe and characterize the charge fluc- 
tuations and correlations of periodically driven single- 
electron emitters on all important time scales. As a spe- 
cific example, we have considered a mesoscopic capaci- 
tor for which we demonstrated that the WTDs contain 
considerable additional information compared to what is 
available in the FCS alone. 
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